Tunneling resonances in quantum dots: Coulomb interaction modifies the width 
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Single-electron tunneling through a zero-dimensional state in an asymmetric double-barrier 
resonant-tunneling structure is studied. The broadening of steps in the I-V characteristics is found 
to strongly depend on the polarity of the applied bias voltage. Based on a qualitative picture for 
the finite-life-time broadening of the quantum dot states and a quantitative comparison of the ex- 
perimental data with a non-equilibrium transport theory, we identify this polarity dependence as a 
clear signature of Coulomb interaction. 
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Single-electron tunneling through zero-dimensional 
states has been observed in a wide variety of systems, 
including metallic islands, lateral quantum dots in gated 
semiconductor devices, vertical dots in double-barrier 
resonant tunneling structures, and molecular systems 
such as carbon nanotubes 0. The I-V characteristics 
has the shape of a staircase, in which each step is as- 
sociated with the opening of a new transport channel 
through the system. Manyfeatures observed in the trans- 
port measurements |M 1^ Q IM IS can be explained either 
within a single-particle picture for non-interacting elec- 
trons or by the orthodox theory of sequential tunneling 
0, H, valid for weak dot-lead tunnel coupling. This is 
no longer the case for interaction effects on transport be- 
yond the weak-tunneling limit. A famous example is the 
zero-bias anomaly of Kondo-assisted tunneling [l^ . 

In this work, we report on a new clear signature 
of Coulomb interaction beyond weak tunneling, that is 
achieved under much less stringent experimental condi- 
tions than required for the Kondo effect to occur. This 
signature is contained in the width of the first step of the 
I-V characteristics. We observe that the width strongly 
depends on the polarity of the applied bias voltage. This 
behavior can neither be explained within a single-particle 
picture nor by sequential-tunneling theory Due 
to Coulomb interaction, the finite-life-time broadening of 
the dot levels becomes energy dependent. As a conse- 
quence, the values for the broadening at the two consid- 
ered steps of opposite polarity can differ by up to a factor 
of two for strongly asymmetric coupling strengths of the 
two tunnel barriers. We use the results of a diagram- 
matic real-time transport theory tliat includes the 
above described physics to find reasonable quantitative 
agreement with the experimental data. 

The experiment was performed with a highly asymmet- 
ric double-barrier resonant-tunneling device grown by 
molecular beam epitaxy on an n^-type GaAs substrate. 
An undoped 10 nm wide GaAs quantum well is sand- 
wiched between 5 and 8 nm thick Alo.3Gao.7As tunneling 
barriers separated from highly-doped GaAs contacts (Si- 
doped with ngi = 4 X IQ-'^^cm"'^) by 7 nm thick undoped 
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FIG. 1: (a) Schematic energy diagram of the asymmet- 
ric double-barrier device under finite bias, (b) First cur- 
rent step in negative bias direction, (c) Comparison of the 
full width half maximum (FWHM)-value of the differential- 
conductance peaks for both bias polarities at base tempera- 
ture (T=20 mK). 



GaAs spacer layers. The sample was fabricated as a pil- 
lar of 2 yum diameter. Two-terminal dc-measurements 
of the I-V characteristics were performed in a dilution 
refrigerator at temperatures between 20 mK and 1 K. 
The studied GaAs quantum well embedded between two 
AlGaAs barriers can be viewed as a two-dimensional sys- 
tem with the edges and residual impurities confining the 
lateral electron motion and thus forming dots. Tunneling 
through the energetically lowest state of the dot, at the 
energy Eq = 33 meV [ij] and with a lateral extent of 10 
nm, produces the lowest resonance peak in the differen- 
tial conductance G = dl/dV. 

The system considered is sketched in Fig. ^a) . Elec- 
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trons tunnel from the heavily-doped emitter through the 
spin-degenerate quantum-dot level embedded inside the 
quantum well. The level comes to resonance with the 
emitter's electro-chemical potential at a finite bias volt- 
age, indicated by a step in the current, as shown in 

Fig.inb). 

Coulomb interaction prevents double occupation of the 
level for the considered range of bias voltage. This gives 
rise to charging effects that are differently important for 
the two polarities. The bottleneck of transport is pro- 
vided by the thicker tunnel barrier. For V < 0, the 
"charging direction" , as sketched in Fig. ^ (a) , the dot 
is predominantly singly occupied, and, therefore, the two 
spin channels effectively block each other. This reduces 
the current-step height to one half of the value for the op- 
posite polarity, V > 0, the "non-charging direction'' for 
which the dot is predominantly empty i^, iSl, S ll^ ■ 
This is well understood within the sequential-tunneling 
picture. 

The focus of this paper, however, is the finite broad- 
ening of the step edge. This broadening can be mea- 
sured as the full width at half maximum (FWHM)-value 
of the differential conductance (dl/dV) peaks associated 
with the current steps. Figure^c) shows our experimen- 
tal low-temperature differential conductance peaks for 
both bias polarities. As a result, the width AV- = 152 
/iV for negative bias has roughly twice the value than 
that of AV+ = 71 /iV. In previous experiments the res- 
onance width in asymmetric structures has been already 
studied and a polarity-dependence has been observed, 
but their explanations referred to polarity-dependent 
leverage- factors 3] or internal saturation processes [3. 

The width of the current step edge at low tempera- 
ture reflects the finite-life-time broadening of the zero- 
dimensional state due to tunneling in and out of the dot. 
The golden-rule rates for an electron tunneling in and 
out of the dot are given by r+(w) ~ J2r=LR^rfr{t^) 
and r~(w) = X]r=LRrr[l — fr{i^)], respectively. The 
tunnel-coupling strength is characterized by the constant 
Tr = 2i:Vr\tr\'^ , whcrc is the density of states in lead r, 
and tr is the tunneling amplitude, and /r(w) = f{uj — fir) 
is the Fermi function of lead r with electro-chemical po- 
tential fir- We note that the presence of strong Coulomb 
interaction introduces an asymmetry between tunneling- 
in and tunneling-out processes: while for an empty dot 
there are two possibilities to choose the spin state of 
the incoming electron, the spin state of an electron leav- 
ing the dot is fixed. It is, therefore, the combination 
2T+{uj) + T-{u), i.e.. 



^F, [1 + /.H] 



(1) 



that determines the finite-life-time broadening. This ex- 
pression is energy dependent, i.e., it depends on the rela- 
tive position of the relevant transport channels to the 
Fermi energy of the leads. At low temperature, the 
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FIG. 2: Sketch of our device in non-charging direction. The 
electronic structure of the dot is given by the spectral den- 
sity p{uj), which is shifted and broadened with respect to the 
bare level eo -I- (1 — a+)eV (solid line) as reflected in real and 
imaginary part of the self-energy E, respectively. 



broadening probed by the electrons within the transport 
window set by the Fermi energies of emitter and collector 
is either 2Fl + Fr or Fl + 2Fr , depending on whether the 
left lead serves as emitter or collector. This qualitatively 
explains the polarity-dependent step width. It is an inter- 
action effect since in the absence of Coulomb interaction 
the asymmetry between tunneling-in and tunneling-out 
processes is lifted by processes involving double occu- 
pancy of the dot, and the finite- life-time broadening is 
energy independent, given by Fl -I- Fr. 

For a more quantitative analysis we employ a dia- 
grammatic non-equilibrium transport theory In par- 
ticular, we make use of the results for the current ob- 
tained within the so-called resonant-tunneling approxi- 
mation for transport through a zero-dimensional state in 
presence of strong Coulomb interaction such that double 
occupancy is prohibited. The technique, the approxima- 
tion scheme, and the steps of the calculation are pre- 
sented in Ref. The result for the current is 
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[uj -e-Rc ^iuj)f + [Im ' 



(2) 



where e is the dot level energy [with the appropriate 
incorporation of a bias voltage influencing level energy 
e = Eo + (1 — a+)eV and Fermi-energies of the leads: 
/iR — const. < eo, /zl = Ml + ^'^^ ^he non-charging 
direction] and the self energy (see Fig. ^ 



Evaluating the integral leads to Rel](aj) 
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and 



lu,^{0j) = -Y^^[l + fr{ 



(4) 



The real part is weakly dependent on a high-energy 
cutoff Ec given by the smaller of the charging energy 
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FIG. 3: Experimental (solid lines) and theoretical (dashed 
lines) differential conductance for the (a) charging and (b) 
non-charging direction at base temperature (T=20 mK). 



for double occupancy or the band width of the leads. 
In the imaginary part, we recover the structure of the 
finite-life-time broadening as postulated in the qualita- 
tive discussion above. Deep in the Kondo-regime the 
approximation above is no longer valid. The two spin 
channels become independent from each other and 
ImE(c^) = -E,r./2. 

We remark that for the step height the well-known 
(sequential-tunneling) result A/-|_ = (2e/?i)rLrR/(2rL + 
Fr) for the non-charging direction [and L <-> R for the 
charging direction] is reproduced. 

The focus of this paper, however, is on the width of the 
step edge. At low temperature, Eq. Q simplifies to 
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with e = eo + {l — a^)eV and q:+ being the leverage factor 
for positive bias [and a- = 1 — a+ being the leverage- 
factor for negative bias], denoting the voltage drop over 
the left barrier. Neglecting the real part of the self-energy 
for the moment, we find that the differential conductance 
as a function of is a Lorentzian with a FWHM of: 

a+eAy+ = 2Fl + Fr, and a_eA\/_ = Fl + 2Fr . 

For a+ « Q!_ and strongly asymmetric tunnel-coupling 
strengths, Fl < Fr, we get the relation AV~ w 2AV+. 

At high temperature we find for the FWHM eAV± = 
3.525A;BT/a± + eAV±., i.e. the temperature broadening 
of the Fermi-function. As predicted from sequential- 
tunneling theory the width increases linearly with 
temperature. The constant term AV^ is of the order 
of F, and, in general, also polarity dependent [l^ . 

For a detailed comparison between theory and experi- 
ment we need first to determine the system parameters. 



The factor a determining the bare level shift with bias 
voltage is gained from the linear high-temperature de- 
pendence of AV±{T) as a+ = 0.53 and a_ = 0.5 so 
that a-f + Q!_ w 1. The coupling constants Fl/r could, 
in principle, both be determined from the current steps 
AI± . For strong asymmetry, however, as is the case here, 
the maximum current is limited entirely by the bottle- 
neck of the smaller coupling Fl. As a consequence, the 
step heights only fix Fl ^ A/+ft/(2e) = 2AI-h/{2e) = 
0.64 ^eV. 

Considering the results derived above for zero temper- 
ature FWHM, we note that this width is essentially the 
sum of different couplings and, hence, is dominated by 
the larger coupling Fr. Therefore we can find Fr = 
{\-a+)eAV+ w {\-a-)AV- /2 . In fact, we gain better 
accuracy by fitting the full dl /dV-^eak for charging po- 
larity and lowest temperature to pinpoint Fr — 40 /xeV. 
The high-energy cut-off Ec = 30 meV is given by the 
bandwidth, i.e. the value of the Fermi energy, being of 
the same magnitude as the charging energy. 

Figure |21 shows the comparison between experimen- 
tal data and theoretical calculations for the differential 
conductance. For the charging direction (negative-bias). 
Fig. |3fa), we find a good agreement between experiment 
and theory for both the resonance width and amplitude. 
In the non-charging direction (positive bias). Fig. Efb), 
the peak width is reduced by about a factor of two for 
both experiment and theory. The experimental data 
show some extra features. First there is an oscillatory 
fine structure on the positive-voltage side of the reso- 
nance attributed to the fluctuations of the local density 
of states of the emitter, see, e.g., Ref. 0. Moreover, we 
observe an enhanced resonance amplitude as compared to 
the theoretical calculation. This effect may be related to 
an additional many-body phenomenon at the Fermi-edge 



The temperature dependence of the resonance width 
in the broad range between 20 mK and 400 mK is shown 
in Fig. ^a). The polarity dependence of the width is 
clearly visible. At low temperature experimental and 
theoretical data nicely match. For temperatures above 
200 mK the width increases linearly with ea±AV± — 
2>.b2bk-QT + C(F). This experimental result is again ad- 
equately reproduced by theory. 

Further support for our explanation of the polar- 
ity dependence as an interaction effect is given by the 
magnetic-field dependence. The data, shown in Fig.^b), 
reveal that polarity dependence is strongly reduced by a 
Zeeman-splitting (on a scale of « 1 T, which is of the or- 
der of the magnitude of the coupling) . In this case, only 
one spin state contributes to transport, and the system is 
equivalent to a noninteracting one, for which theory pre- 
dicts Iml](a;) — — Fr/2 and, thus, a low-temperature 
width a±eAV± = Fl + Fr independent of the polarity. 
This is in accordance with the trend seen in the experi- 
mental data. 
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FIG. 4: (a) Temperature dependence of the resonance width 
AV for both polarities, open circles V > 0, filled squares V < 
0, continuous lines theoretical simulations, (b) Saturation 
width as function of B-field for both polarities. 



The discrepancies remaining between experimental 
data and theoretical simulation are consistent with the 
expected range of accuracy. The main source of ex- 
perimental errors are systematic fluctuations in the cur- 
rent, in particular the local DOS fluctuations prominent 
in non-charging direction which become even more pro- 
nounced at high temperatures. This artefact particularly 
limits the precision to which we can determine the cur- 
rent step height A/+. 

We study the width of the current step of single- 
electron tunneling through a zero-dimensional state in 
a double-barrier resonant tunneling structure and flnd 
that it is strongly polarity dependent. This is interpreted 
as a clear signature of Coulomb interaction. The lat- 
ter introduces an asymmetry between possible tunneling- 
in and tunneling-out processes, which gives rise to an 
energy-dependent finite-life-time broadening. The bias- 
polarity dependent step width, observed in our exper- 
iment, where strengths of the tunnel coupling between 
dot and the two leads were highly asymmetric, can be 
simulated with reasonable agreement by applying a di- 
agrammatic non-equilibrium transport theory for inter- 
acting quantum dots. 
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sistance was granted by BMBF and by DFG via SFB491 
and GRK726. 
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